CRITICAL ANGLE IN THE PARTON CASCADE 



Wolfgang Ochs 
Max- Planck- Institut fur Physik, Werner Heisenberg Institut, 
Fdhringer Ring 6, D 80805 Miinchen, Germany 

Jacek Wosiek 
Institute of Physics, Jagellonian University, 
PL 30-059 Cracow, Reymonta 4, Poland 



Abstract 



The angular correlation function of partons in a jet as derived from perturbative QCD is 
nonanalytic at a critical angle which separates a multiparticle regime with scaling properties 
from a regime with few particles near the hadronisation scale. Some phenomenological 
consequences are discussed. 



MPI-Ph/95-98 
TPJU-17/95 
November, 1995 
hep-ph/9511457 



1 



I. INTRODUCTION 



Recent studies of angular correlations of partons inside jets within pertubative QCD have 
demonstrated the emergence of scaling properties at high energies and at sufficiently large 
relative angles [|l|- ||. In this paper we focus on the complementary kinematic region of small 
relative angles within the same theoretical approach of perturbative QCD. The correlation 
function becomes nonanalytic in the high energy limit with a discontinuous second derivative 
for a critical angle d c which separates two regimes of quite different characteristics. In the 
small angle regime the correlation function is determined by the one-particle energy spectrum 
of the parent parton from which the detected partons emerge. In the theory with running 
a s the correlation function in this new regime depends on the nonperturbative transverse 
momentum cutoff Qq in the cascade whereas in the large angle regime there was no cutoff 
dependence. This opens the possibility to test the hypothesis of "Local Parton Hadron 
Duality" (LPHD) || near the boundary Q . 

The occurrence of a critical angle d c in the angular correlation has been noted already 
in two previous papers. In our first study on angular correlations we found (for fixed a s ) 
a relation between 2-particle correlations, and the 1-particle spectrum for the critical angle 

^(™T (id 

p\qJ 1 

and also for correlations in a pair of points symmetric in certain variables around this critical 
angle (G and P are the opening angle and primary momenta of the jet respectively). 

In the analysis by Dokshitzer, Marchesini and Oriani |l| the same angle ( |1 . 1|) was found 
as limiting the validity of a power law in the relative angle from below (again for fixed a s ). 
Moreover they suggested the existence of a correlation function below d c in a nonperturbative 
model: the rate for the particles with relative angle $12 is given in terms of the momentum 
distribution of the parents. 

We show in the following that perturbative QCD in its regime of applicability, ie. for 
transverse momenta k T > Q , provides the correlation function not only above but also 
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below i? c . It can be expressed indeed by the momentum distribution of the parent parton in 
some analogy to the above effective hadronic model Furthermore we extend the analysis 
to the interesting case of running a s . Again a critical angle can be identified. At high 
energies it behaves as 

^A^qVA)^ (12) 

which depends more weakly on the jet energy than the angle ( |1.1| ) for constant a s (A is 
the QCD scale and A = ln(Qo/A)). The special role of the critical angle as singular point 
becomes obvious from the fact that the correlation function depends on the cutoff Qo only 
below but not above it. 



II. EMERGENCE OF THE CRITICAL ANGLE IN PERTURB ATIVE QCD 

We consider for definitness the correlation density p^ 2 \$i2, P, ©) in the relative angle $12 
of two partons and comment later on multiparton correlations. In the double logarithmic 
approximation (DLA) the correlation function obeys the integral equation H 

P w($ 12 ,p,e) = dv\# 12 ,p,e)+ f P ^ / e ^nm^^K^). (2.1) 

JQo/012 A J0 12 W 

where the anomalous dimension 70 for the evolution of the multiplicity is given in terms 
of the running coupling 7^ = 6a s /vr = (3 2 / In (k T /Q ) with 1 = 12JV C /(11JV C - 2N f ). The 
inhomogenous term is constructed from the product of 1-particle angular distributions 
and found as S 2 \x) = p^(x)n(x), x = Pfiril ^ and with the multiplicity n in the cone of 
half-angle $12 around the primary parton axis. In exponential accuracy S- 2 ^ ~ n 2 . 
The leading contribution to the solution of Eq. fl2.1|) was derived in 00 



Present R{P/K) denotes R{P/K) + 5(1 - P/K) of Ref. §. 
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The resolvent R^ } ^ } ^a 12 ) is the momentum distribution of the parent partons K in 
the jet (P, 0) with the condition, that their virtuality is bounded from below by Kd^ and 
not by the elementary cut-off Qq. This is equivalent to limiting from below their emission 
angle by #12. The representation ( |2.2| ), which was derived as the solution of the differential 
equation (|2.1|) is identical with the jet calculus result used in Ref. PJ. 



In order to show the existence of the nonanalyticity of the p^'^x-z) given by Eq. ( |2.2j ) at 
some small value of the relative angle $12 = consider the dependence of the integrand 
of Eq. (|2.2j ) on the parent parton momentum K. The interplay between the decrease of the 
parent density R and the increase of the children densities with K produces a sharp 
maximum which is however ■i9 12 dependent. The important element of running a s analysis 
is that the position of this maximum is independent of Qo, since both distributions do not 
depend on Q Q f\ Now, for small relative angles the above maximum shifts below the lower 
limit Qo/^12 of the momentum integration. At that point the result changes nonanalytically 
(albeit very gently - analogously to the second order phase transitions in statistical systems). 
Namely, for bigger than p^ is given ( to the exponential accuracy) by the Qo 
independent value of the integrand at the maximum. For $12 below d c however, the result 
is given by the value of the integrand at the lower bound and obviously depends on Q . 
Moreover, since <i 2 (l) = 1 in our exponential approximation, the density of pairs below d c 
is given entirely by the density of parent partons at the minimal momentum K = Qo/'&u- 

It is easy to verify explicitly the above considerations for constant a s . In the logarithmic 
variables which are natural for this problem 

F = ln(P6/Q ), u = hi(Ptf 12 /Qo), v = \n(Q/$ 12 ), (2.3) 
£ = ]n(P/K) =ln(l/x), (2.4) 

the solution (p. 21) reads 



2 Except of the overall normalization factor which does not influence the shape of the momentum 
dependence, hence is not relevant for the argument. 

4 



p^(u,Y)= [ U d^,v)d 2 (u-0- (2.5) 
Jo 



Using the known momentum distribution of parent partons p*- 1 -*^, Y) = exp (27oy £(F — £)) 
and replacing the vituality cut-off Qo — > K$i2, i-e. F — > t> + £, gives 

p (2) («,F) = jT exp (70(2^ + n(u (2-6) 

where we have also substituted the asymptopic form of the inhomogenous term G? n (w — £) = 
exp (ri7o("U — £)). The integer n, which for the two-body correlations is just 2, will be kept 
arbitrary in the following discussion. This allows for the direct generalization to the higher 
order correlations and multiplicity moments. The maximum of the integrand in fl2.6| ) is 
located at 

Cmax • (2 .7) 

V 



This remains in the integration range only if 

n 2 + l 



?<^T7-r, (2-8) 



where we have introduced the scaling variable 

v _ ln(8/0 12 ) 
u + v ln(P6/Q )' ' 

In this case the result for the density reads 

p {2 \u,Y) ~ exp f 7o (n« + -)Y (2.10) 



and we obtain the power behaviour in the angle p^ ~ i?^ 70 ^- If £,max > u the integral i 



is 



saturated by the maximal value of the integrand at the boundary and we obtain 

(u, Y) ~ R(u, v) ~ exp (2 7ov W) , e > , (2.11) 

v ' n z + 1 

i.e. the two-parton distribution is given in terms of the single parton momentum distribution. 
One can summarize our results in the "e - scaling" form 



5 



P {2) (e,Y) ~ exp( 7o Y£(e)), (2.12) 
w(e) = n - (n - ±)e, e < e c , (2.13) 



Q(e) = 2y/e(l-e), e>e c . (2.14) 

where e c = n 2 /(n 2 + 1). Note that the solution ( j2.14j ) is n— independent above e c . The 
critical angle in Eq. ( |1 . 1| ) quoted in the Introduction corresponds to e c = 4/5 for n — 2. It 
is always bigger than the minimal relative angle "dram — Qo/P and the maximal transverse 
momentum associated with $ r grows like (P/Qo) 1 / 5 . The same result can be obtained from 
the exact solution for p^ in terms of Bessel functions |2| which is valid also at finite energies. 
The above solutions are displayed in Fig.l for different n. The linear solutions for e < 



e c are tangent to the ellipse ( 2.14 ); at the critical point e c the functions Coil) and their 
derivatives are continuous but not their second derivative. 

The same asymptotic result ( 2.12 ) for p^ is also obtained from the nonlinear, first order 



differential equation 

(u - eu>')(u> + (1 - e)u') = 1 (2.15) 

which corresponds to Eq. fl2.1|) ||. One finds the linear solution, Eq. Q2.13|) , for any boundary 
value oj(0) = n > and the ellipse solution, Eq.( [2.14| ), for the boundary value cD(l) = 0. 
Eq. (|2.15|) can be solved explicitly for Co' with two roots in general; for a point on the ellipse 
there is only one root corresponding to a bifurcation point where the two solutions with 
proper boundary conditions at e = and e—1 meet. 

The simpler case of fixed a s or 70 considered above is instructive as it allows for the 
exact solution 0, in this case the dependence on the QCD scale A in Eq. ([2.1| ) drops out 
and the only reference scale is provided by the cut-off Q$. The constant a s solution is also 
obtained from the running a s result when both (3 2 and A = InQo/A — > 00 with their ratio 
/9 2 /A = 7q fixed. This is equivalent to the famous limit Nf — ► llN c /2 with A — > tuned in 
such a way that the anomalous dimension 70 remains fixed. In other words A behaves as 

A = ^«p(- ^ Uo, (2.16) 



in this limit. 



III. RUNNING a s AND DOUBLE SCALING 

We employ the Mellin representation of the inclusive momentum distribution for running 
a s , see Appendix B of Ref. [0] for the details. With the appropriate asymptotic forms of 
the energy moments one obtains 

p^($ 12 ,Q,P) « / exp(sv)F(s,L)exvw n (s,l)^, (3.1) 



where 



2/3 2 sfx . 



(3.2) 



s y 2(3 1 

expw n (s,l) = f l x da expv n (s,a), (3.3) 

— In • 

s y 2p 



v n (s, a) = - %g-(s, a) —In (—g + (s, a)) + 2nf3</a, 



g±(s,x) = ^/s 2 + A/3 2 /x±s, (3.4) 

and all prefactors have been neglected. The logarithmic variables are now defined with A as 
the reference scale 

L = ln(P0/A), Z = ln(Ptf 12 /A), a = In (K# 12 /A), A = ln(Q /A), (3.5) 



and v as in Eq. fl2.4|) . The saddle point of the a integration is located at 

, /3 2 (l-n 2 ) 2 

= 7w ■ (3 ' 6) 

However the maximum of the a integrand is not always within the bounds of the real a 
integration. In particular for 

s = s Q = — j= , (3.7) 

this maximum coincides with the lower cut-off at a = A. Second characteristic point in 
the complex s plane is the saddle point s x of the s integration. Depending on the external 



parameters (e.g. $12) s x can be bigger or smaller than s thus determining the leading 
contribution to the integrand [] . The condition that the two points coincide determines the 
critical angle $ c . In summary: for all $12 the density of pairs can be written in the form 

( e , L, p) ~ exp (2pVZ(u n (e, p) - 2p)) , (3.8) 

with 

v _ In (g/gu) _ /777 ^ 

e -Z-ln(P9/A)' ^-VV^- (3-9) 

For i?i2 > $ c both integrations are done by the saddle point approximation with the result, 

z = sVL/{3, 

u n {e) =j(z(e))+ez, (3.10) 



j(z) = -(y/z* + 4-z), (3.11) 

and z n (e) is the solution of the following algebraic equation 

7 2 (z) -21117(2) - ez 2 = n 2 - 2 Inn. (3.12) 

For i?x2 < i?o one uses the value of the a integrand at the lower boundary a = X as the 
estimate for the a integration which results in the different saddle point of the s integration. 
The final result reads in this case 

u n (e, p) = 7(2) + ez - p(>y(zp) - n), (3.13) 

and z n (e,p) is determined by the different saddle point condition 

7 2 (^) - 2 In 7(2) - ez 2 = ^ 2 (zp) - 2 In 7(3,0). (3.14) 

The critical point 

z Q = s e VL/{3 = -^±, (3.15) 
np 



3 The value of s x is real. 
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coincides with the saddle point at $12 = $ c or equivalently at e = e c . This fact can be used 
to determine e c from e.g. Eqs.( |3.l2J ) and (|3.15| ). The consistency check of such a procedure 
is the equivalence of both equations ( |3.12| ) and ( p. 14 ) at z = z & . This is readily verified 



since 7(z p) = n. Our result for e c reads 

e c (p) = \ (l 2 {z & )-n 2 - ln( 7 2 (^ )/n 2 )) . (3.16) 



Zq 

or explicitly for n = 2 (see Fig. 2) 



e c ( P ) = \qi+ (f ) 2 + I) 2 - (f ) 2 - f In ( ^tyl + (f ) 2 + l)/2 ] . (3.17) 

Figure 3 shows the u functions obtained from Eqs. (|3.10| - |3.14j) for few values of p . As 
expected from the general discussion in Sect. 2, e c and u n (e, p) (for e > e c ) depend on Q . 
This is in contrast to the constant a s case, where the only scale was provided by Qo itself. 
For the running a s , the second scale A allows to construct two dimensionless ratios relevant 
for this problemfj, and, as a consequence, the simple e scaling is violated for e > e c ($12 < $ c )- 
On the other hand, u n (e, p) is independent of n for e > e c (n) as in the fixed a s case. 

It follows from Eq. flOp that the new solution, Eq. (|3.13j ), has the correct behaviour 
at the kinematical boundary $12 = Qq/P. One verifies readily that Eq. (|3.14|) has a real 



solution only for e < e rnax = 1 — p 2 which corresponds to the above kinematic limit for ^ 12 . 
Around that point uj has the following expansion 



u(e,p) ~ 2p - 2\n(p)s/(e max - e) + A 3 ^J (e max - e) + ... (3.18) 

with 

a = 2p 2 In 3 / 2 (1/p 2 ) - 2p 2 In (1/p 2 ) - (1 - p 2 )(l + 2p 2 ) 

3 2p 2 ln 3/2 (l/p 2 ) ' 1 ' ' 

This is to be contrasted with the asymptotic behaviour of the infinite energy limit u(e) at 

e ~ 1 obtained from Eqs.( [3.10| - |3.12| ). 



4 In the DLA the limit A — > 00 would lead to a divergence of the multiplicity and is therefore not 
allowed 
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W ( e )~y(l- C )]n-U (3.20) 

The finite value of u(e, p) = 2p at the boundary ensures our proper high energy normalization 
of the density p^ ~ 1 at #12 = Qo/P. The coefficient ln(l/p 2 ) of the threshold behaviour 
turns into the additional logarithmic singularity of the asymptotic solution (|3.20| ) around 
e ~ 1. 

The constant a s results of the previous Section can be recovered from Eq. (|3.16|) in the 
limit A — > 00 or p — > 1. One finds 

e c ^(l-p 2 ), p-1, (3.21) 



which agrees with ( |2.14 ) since the scaling variable for the constant a s , e = e/(l — p 



2\ 



On the other hand in the high energy limit, L — > 00, Q — fixed, p — > and Eq.( p,17| ) 
gives 

4p 2 , , 1 N , 

ln(-), p^O, (3.22) 



and the critical angle 
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c = pexp(L(l-e c (p))), (3.23) 



may become small raising the question of the applicability of perturbative QCD in the new 
regime ($12 < i? c )- However it is easy to show that the transverse momentum Pd c is always 
bigger than the perturbative cut-off Q which can be written as 

Qo = Aexp (Lp 2 ). (3.24) 

Since 

e c (p) < e maa: = 1 - p 2 (3.25) 

is always in the physical region ( d c > d™ n — Qo/P ) , perturbative QCD applies also for 
Qo/P < $12 < The above argument is valid only for p > 0. Exactly at p = (infinite 
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energy, at fixed Q$) the inequality (|3.25|) is saturated and using the asymptotic form (|3.22|) 



we obtain result ( |1.2j ) quoted in the Introduction. It implies the transverse momentum 
associated with i9 c grows as a power of the logarithm of P and is therefore larger than Qq. 

At the critical angle the u> function, hence also p^ 2 \ are known analytically since the 
solution of the Eqs.( ^.l2] , ^.14D , i.e. z & , is known explicitly, c.f. Eq.( ^.15[ ). In particular in 



the high energy limit, Eq.( |3.22 ), we obtain 

lnfpe/A; 



M t „ P ,e ) = ( , 26 , 



A 

i.e. the density of pairs with the critical opening angle is growing only with a power of the 
logarithm of the energy. 

The transition between the high energy and constant a s limits can be conveniently 
presented in slightly different variables. In Fig.4 we show the function Wjj(e, p) 

cu R (e,p) = (u(i(l-p 2 ),p)-2p)/(l-p), (3.27) 

as the function of e, in the normalized interval < e < 1, and for several values of p. The 
constant a s regime (c.f. Fig.l ) is already reached for p ~ .8 

The limiting behaviour, Eq.( |3.22|) implies that the whole region $ 12 < $ c shrinks to zero 
(e = 1) at infinite energy and constant Qq. Therefore the existence of the non-scaling part 
u(e, p) can be regarded as the non-leading effect. However there is an important distinction 
between this type of finite energy corrections and the behaviour for $12 > $ c , namely below 
d c the leading contribution is Qq dependent, while above d c it is not. This follows from the 
general argument presented in Section 2. Second point is that precisely these corrections 
lead to the finite value of e c in the constant a s limit. In other words, the high energy and 
constant a s limits are not interchangeable at fixed Qq, and the region -i9 12 < $ c is responsible 
for the difference. With that point in mind, one can define the double scaling limit: 

L — > 00, e, p — fixed. (3.28) 

In that limit u(e, p) is the true asymptotic scaling function of two scaling variables, for all 
< e, p < 1. It is singular at e = e c (p) and becomes independent of p for e < e c (p). The 
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double scaling limit ( |3.28| ) requires Qq to grow with energy ( L ) according to Eq. ( |3.24| ). 



This is not quite along the standard assumption of constant Qq at high energies. However 
this limit allows us to make precise prediction in the L —>■ oo limit for the two regions above 
and below e c which are expected to persist for finite energies. 



More importantly, comparison of Eqs. ( |3.23| ) and ( |3.24| ) shows that the necessary condi- 
tion for applicability of the perturbative QCD, is now safely satisfied since the Pd c grows 
exponentially faster than Q due to (|3.25|) . Therefore even the softer parents are described 
by the perturbative approach in this limit. The singularity nature of the u function at e c (p) 
follows readily from Eqs. (|3.10| - |3.14|) : u)(e, p) is continuous together with its e derivative but 



the second derivative is not. 

By changing p within < p < 1 one can move the bifurcation point e c through the whole 
range of < e < 1. Since at e c the uo function is known analytically and, moreover, the 
solution of the algebraic equations ( |3.12| , |3~7L4|) is also known (c.f. Eq. (|3.15|) ) we can provide 



the explicit analytic form for the scaling function in the parametric representation for all e 



1 / /~. "77 . \ o . o ^' 
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e(t) = -(VT+& + l) 2 - t 2 - - In (^(VT+^ + 1) J (3.29) 

u(t) = --e(t) + -Ul +t 2 + 1), < t < \. (3.30) 

For < t < 4/3 the asymptotic curve oo(e,p = 0) is produced. For finite non-zero p the 
representation (|3.30|) is valid for 4p/3 < t < 4/3 and describes the p independent branch of 



the scaling function U) n (e) for e < e c . 



IV. CONCLUSIONS AND PHENOMENOLOGICAL APPLICATIONS. 

The angular correlation functions exhibit two kinematic regimes at high energies with 
quite different properties, separated by a critical angle $ c or e c . Whereas for fixed a s this 
critical behaviour occurs for fixed energy independent e c , in the realistic case of running a s 
we find e c — >• 1 for P — > oo. We consider here the preasymptotic region with e c < e max < 1 
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(for Q /A > 1) where besides the leading terms of the order In ■s/L, the terms of all 



orders in \J\/L are kept (A = \n(Q /A), L = ln(P6/A)). Then the boundary conditions 
with finite e max can be fulfilled. The occurence of the critical value e c is quite special to 
angular correlations and, for example, there is no such critical value in the distribution of 
the energy variable £ = Inl/x (the "hump backed plateau"): so, in case of fixed a s , the 



function In p^(£) is asymptotically given by y£,(Y — in the full kinematic region whereas 
the angular correlation function is built up from the two different pieces in Eq.( [2.14|) with 
discontinuous second derivative. 

In the following we want to point out two differences between the two regimes (in the 
theory with running a s ) with their phenomenological consequences. 
1. Dependence of the cutoff Qo- 

In the region e < e c the correlation function is Qo-independent and develops a scaling 
behaviour (e-scaling ||) whereas in the complementary region (e > e c ) this scaling is broken 
by Qo). This leads to two possible effects: 

a) as a function of energy the small e region approaches the scaling limit which depends on the 
parameter A, but not on Q , whereas at large e the correlation function drops towards e max 
which depends in addition on the cutoff Qo with e max — > 1 for P — > oo, so it is nonscaling. 
This opens the possibility to determine independently the parameter A from the small e 
data and the effective cutoff Qq from the large e data. The first results from DELPHI 
on "e-scaling" could actually be interpreted along these lines: data with smaller PQ 
drop down earlier with increasing e. For more quantitative results an analysis with higher 
accuracy seems to be required: on the theoretical side an analysis beyond the DLA would 
be desirable and on the experimental side a better treatment of the jet axis problem, for 
example, by using the Energy- Multiplicity-Multiplicity observable ||||; 

b) another effect is the dependence of the correlations on the species of particles under the 
assumption that the particle mass acts as an effective transverse momentum cutoff and is 
therefore related to Qo PI- The correlation functions for particles of different species (say 
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7T7T, KK, pp) should then be the same for small e and different for large e where for 

heavier particles vanishes for smaller e max . 

2. Dependence of the correlation functions on their order n. 

A convenient measure of correlations of higher order are the multiplicity moments, for ex- 
ample the factorial moments =< n(n — l)...(n — q + 1) > in a limited angular ring 
(dimension D — 1) or cone (D = 2) of half opening 5 at polar angle to the primary 
parton (jet) direction 0-0. These moments (usually one considers the normalised moments 
p(g) = f{i) J < n behave at high energies likef] 

In F (n) ~ 2[3VLu n (e, p). (4.1) 

One can determine the function u n (e, p) from the measurement of F n after appropriate 
rescaling ("c^"). Again for smaller e one should obtain the predicted scaling limit 
uj n (e, 0) whereas at sufficiently large e > e c (n) there is the dependence on p = JX/L but 
no dependence on the order n. This prediction follows as the opening angle S determines 
the parent momentum K > Qo/5 independent of the number of particles inside the cone 
because = 1 at the threshold independently of n. In the same way one can extract 
^(e, p) from p( 2 \ Then all rescaled correlations should yield the same scaling function 
ui(e, p) independent of n for large e 

(4.2) 

Also in this region (large e) the cutoff Qq could possibly be determined as lower limit of K5 
or Kdi2 where K is the momentum of the n— cluster or the pair considered. 

We have discussed here the consequences of perturbative QCD in a preasymptotic region 
at high energies. At the energies realistic today one has to consider also finite energy effects. 
At high energies the width a 2 of the peak in the variable C,/Y of the integrand in Eqs. 

5 Factorial F^ n ' or cumulant moments approach the same asymptotic limit but the nonleading 
corrections appear to be smaller for the HQ. 
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( p.2|J2.5|) decreases like l/\/L and the two regimes become separated at e c . At finite energies 
there is an overlap between these regimes of this order. To illustrate the finite energy 
effects we show in Fig. 5 the rescaled normalized correlation function with the asymptotic 
behaviour [[| 

^ M^) = 2 ^ _ 2VT ^- e) (4 3) 

/in P6/A 



and the exact numerical solution of the integral equation (|2.1|) for at LEP energies 

for different opening angles 0. One can see that the scaling behaviour is approached quickly 
for e < 0.4 whereas for large e the separation of curves is clearly visible. In comparison to 
our analytical calculations for finite X/L (but still L — > oo) the numerical results show in 
the central region a maximum of the correlation of lower height and position in e. It should 
be noted, that the finite d c effects become more pronounced, hence easier to observe, at 
lower CM energies, for example at y/s = 5 and 20 GeV we find e c = 0.65 and e c = 0.75 
corresponding to i? c = 10° and $ c = 2.5° respectively. 

In summary, the angular correlations - different from the case of the well studied energy 
distributions - show two regimes which are dominated either by the scaling behaviour of the 
well developed cascade or the hadronisation effects near the transverse momentum cutoff Q Q 
respectively with specific phenomenological consequences for both regions. 
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APPENDIX A 

In this Appendix we study in detail the analytic structure of the Mellin representation 
of the two-parton density p^ in the constant a s case. This serves as the basis of the gener- 
alization for the running a s case performed Sect. III. The energy moments of the momentum 
distribution 
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f jL n (Y) = p-» ( P k n p^(k,P,e)dk= ( Y e- n yp^(y,Y-y)dy, (4.4) 
JQo/e Jo 



read for constant a 



s 



1 / Y i 

I ( (V^ 2 + 47o +n)exp[-(^n 2 + 4 7o 2 -n)] 

F 



+ (^n 2 + 4 7o 2 - n) exp [--( ^n 2 + 4 7o 2 + n)] ) (4.5) 



The distribution of parent partons required in Eg. ( j2.2| ) can now be reconstructed as 

ds 



R(y,x)= J JL e a n(s,x + v) (4.6) 

Therefore Eq.( |2.2p can be rewritten in the form 

p (2) (e,y) = e^° y(1 - f) J^ dsexp {{eY/2)^s 2 + 7o 2 + s)) (4.7) 
^ y(1_£) dyexp (( 2/ /2)( v /s 2 + 7o 2 - s - 2n 7o )) (4.8) 

where all preexponetial factors have been omitted. The new contour 7 follows from the 
change of variables s — > — s and from the interchange of the s and y integrations. The 
integration over the parent momentum (y = \n(P/K)) is now elementary and one is left 
with 



pW(e, Y) = exp (n 7o y(l - ?)) /, &exp (eY/2)(Js* + 7 „ 2 + s) (4.9) 



exp^((l-e)y/2)(^ + 7g- S -2 n7 o)j ^...^ ■ < 410 > 

where the nonleading contribution from the lower limit have been neglected. The final result 
is determined by the relative positions of the pole and the saddle point of the integrand. 
This analytical structure of is shown in Fig.6. The pole is located at 

1 - n 2 

s& = 7o^—, (4.11) 

while the saddle point position 
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1 — 2e , 
=7o I -- 4.12 

m - e 



depends on e, and the contour 7 is now to the left of these singularities. For e < e c ( Fig. 6a ) 
the saddle point is separated from the contour by the pole, and while deforming the contour 
7 — > 7' one picks up the pole contribution (which is in fact dominating over the saddle 
contribution). This gives the result (|2.12| , |2.13|) . For e > e c however s x < s Q and the pole 



does not contribute, see Fig. 6b. In this case the saddle point saturates the integral and one 



gets the n— independent result (|2.14|) . Analogous structure occurs in the running a s case as 
discussed in Sect. III. 
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FIGURES 




£ 

FIG. 1. Scaling function <D n (e) for fixed a s . For e < e c two solutions corresponding to 
and n = 4 are shown; for e > e c the solution is given by the ellipse for all n. 




P 



FIG. 2. Critical and the maximal values e c (p) and e max (p) as the function of p. 
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FIG. 4. Scaling plot of the function cjfj(e, p) of Eq.( |3.27 ) in the e variable for 
p = 0.001,0.4,0.8,0.999 (from top to bottom); p ~ corresponds to the running a s , high en- 
ergy limit, p rsj 1 to the constant a s limit. 
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FIG. 5. Rescaled normalized correlation function r(e), c.f. Eq. (|4.3| ) for asymptotic energies 
(full lines) and for LEP energies with jet opening angles (a) = 15° and (b) = 60°. The dashed 
curves correspond to our asymptotic predictions of Eq. ( |3.12| , |3.14] ) , the dash-dotted curves to the 
numerical solutions of Eq.([2.1|) (for A = 0.15 GeV, Qo/A = 2). The arrow indicates the position 
of the critical angle 




e<e, 




-c v —' w c 

FIG. 6. Analytical structure of the Mellin representation Eq.( [4,10| ) in the complex s plane; (a) 
e < e c and (b) e > e c . 



21 



